CRM-2512 (1997) 1 
solv-int/9710016 



BISPECTRAL OPERATORS, DUAL 
ISOMONODROMIC DEFORMATIONS 
AND THE RIEMANN-HILBERT DRESSING METHOD 



Q\ ■ J. Harnad 



ON 



Abstract. A comparison is made between bispectral systems and dual isomon- 
i odromic deformation equations. A number of examples are given, showing how 

bispectral systems may be embedded into isomonodromic ones. Sufficiency condi- 
tions are given for the construction of rational solutions of isomonodromic deforma- 
ns) , tion equations through the Riemann-Hilbert problem dressing method, and these are 
£SJ ■ shown, in certain cases, to reduce to bispectral systems. 



o 



> 

■ 1. Bispectral Systems vs. Isomonodromic Deformations 

' We begin by listing a number of results characterizing bispectral systems and 

isomonodromic systems, respectively, with a view to comparison. 

' 1.1. Bispectral Systems. 

(1) In bispectral systems, we have a function ip(x, z) of the two variables (x, z) 
that simultaneously satisfies a pair of eigenvalue equations 



o 



(1.1a) Ltp{x, z) = f(z)ip(x 1 z) 

(1.1a) Aip(x 1 z) = <j){x)^jj{x, z), 



where L(x) and A(z) are differential operators in the indicated variables and 
f(z), 4>{x) are parametric families of eigenvalues depending on the other 
variable. Moreover [DG], for a suitable normalization, L(x) takes the form 

(1.2) L(x) = d: + a r - 2 d r x - 2 + --- 

where the coefficients {aj(x)}j = i... r -2 are rational in x and the eigenvalues 
f(z) depend polynomially on z, and similarly for the operator A(z) and 
eigenvalues (f>(x). 
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(2) For L's belonging to rank-1 bispectral algebras [W], the eigenfunction tp 
may be viewed as the t = (t\ = x, 0, . . . ) value of a Baker- Akhiezer function 
ipw(t, z ) expressed, as usual, in terms of a r-function via the formula 

(1.3a) ^=^^Mt- N ) 



(1.3b) 



T W (t) 
1 1 1 

z' 2? 3^' 



The T-function Tvy is rational, and has the form 



(1.4) Tw^Actix + Y^M-zy- 1 ), 

where (X, Z) is a pair ofnxn matrices representing a point in Wilson's 
complexified, prereduced Calogero-Moser phase space C n . That is, their 
commutator is of the form 

(1.5) [X, Z] = -I + fg T 

for a pair of vectors /, g G C" . The t-dependence of tw is thus given by 
flows generated by the commuting family of Hamiltonians 

(1.6) K := -tr(-Z) n , 

which are obtained by applying the bispectral involution (see below) to the 
Calogero-Moser hierarchy. The Baker function tpw at t = (x, 0, • • • ) may 
be explicitly expressed in this case by the formula 

(1.8) tp w (x, z) = e xz det(I r - (X + x\y x (Z + ziy 1 ). 

(3) Solutions can be constructed from suitably defined "vacuum" solutions 
through the application of Darboux transformations [BHY1, BHY2,K, K]. 

(4) The Baker function tpw and the r-function tw for rank-1 bispectral al- 
gebras correspond to an element W of Wilson's adelic Grassmannian Gr ad 
[W]. The existence of a bispectral pair (L, A) follows from the fact that Gr ad 
admits a bispectral involution b : Gr ad — > Gr ad defined by the commuting 
diagram of invertible maps 



1 1 r 1 1 r 



(1.10) w 



w 



Gr ad j Gr ad 

where the map W : U„C„ — > Gr ad is defined by formula (1.8) and bw ■ C n — > C„ 

by 



(l.ll) 



b w :(X,Z)^(Z\X t ). 
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1.2. Isomonodromic Deformations. 

(1) For these systems, we have an invertible rxr matrix-valued function ^(x, z) 
that simultaneously satisfies a pair of equations 

(1.12a) = U(x,z)^{x,z) 

(1.12b) ^^l =V (x,zMx,z), 

where, in general, the rxr matrix-valued functions U(x, z) and V(x, z) are 
rational in z. The consistency condition 



(1.13) 



dz ^' dx ^ 







implies the invariance of the (generalized) monodromy data of the rational 
covariant derivative operator 

(1.14) V z :=^--V{x,z) 

under the deformations resulting from varying x. In general, there may 
be many deformation parameters (t\ = x,t2,- • ■), each with an associated 
infinitesimal deformation operator 



(1.15) V j :=1--U j (t 1 ,t2,--- ,z) j = 1,2,... 



j 

commuting with T> z and amongst themselves, each determining a 1-parameter 
family of isomonodromic deformations. The structure of the corresponding 
deformation equations was analyzed in [JMU, JM]. 

(2) Isomonodromic systems of type (1.12a,b) have a Hamiltonian structure 
[JMU, HI, HTW] and are generated by certain spectral invariants of the 
matrix V . The logarithmic differential of the corresponding r-function is 
given by 

(1.16) rf(lnr) = ^Hjdtj, 

3 

where {Hj}j = i 2, -- are the spectral invariant Hamiltonians generating the 
equations (1.12a,b). 

(3) Solutions to eqs. (1.12a,b) can be constructed from "vacuum" solutions 
by application of the "dressing method" [NZMP, HI] , based on the matrix 
Ricmann-Hilbcrt problem. For suitably chosen vaccua and loop group el- 
ements, this gives rise to solutions having a rational dependence upon the 
deformation parameter. 

(4) The isomonodromic r-functions correspond to certain special elements W £ 
Gr of the general Segal-Wilson-Sato Grassmannian. Moreover, at least 
for the case where the element V(x, z) has n simple poles at finite z and 
tends to a nonsingular finite value X at z = oo, there exists an equivalent 
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representation of the underlying dynamical equations as deformations of a 
second dual isomonodromic family of differential operators [HI] 

(1-17) V x :=^-V(x,X), 

with V(x, A) tending to a nonsingular finite limit Z T at A = oo. The 
corresponding dual matrix Baker function $f(x, A) is a nonsingular n x n 
matrix satisfying a similar system of equations 

(1.18a) A) = t/(x, A)$(x, A) 

(1.18b) =V(x,\)*(x,\), 

where the n x n matrices U(x, A) and V(x, A) are rational in the second 
spectral parameter A. The dual system (1.18a,b) can be related to the 
original one (1.12a,b) by introducing an auxiliary symplectic vector space 
consisting of canonically conjugate pairs (F.G) ofnxr rectangular matrices 
and applying a Hamiltonian quotienting procedure. The resulting matrices 
V and V in the definition of the operators V z and T>\ respectively are given 
by the formulae 

(1.19a) V{x,z)=X + F{Z-zl r )- 1 G T 

(1.19b) V(x, z) = Z T + F T (X T - AI„) _1 G, 

and the corresponding spectral curves are birationally equivalent (with the 
roles of the loop parameter A and eigenvalue z interchanged). These are 
therefore related by the involutive duality map 



(1.20) (X, Z, F, G, z, A) i — ► (Z 1 ,X,F,G,X, z). 

1.3. A simple example. To further indicate the close correspondence between 
bispectral and isomonodromic systems, consider the following example of rank-1 
bispectral systems (cf. [H2]). In the above notation, choose the matrices X and Z 
to be 

(1.21) X 

The resulting bispectral Baker function is then 

(1.22) i>(x,a) = e xz (l-— ' 











s 


i ■ 


z= I 





which satisfies the eigenvalue equations 

^ /a 3 6 5 12\ , 3 ; 
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In this case, we have two other single-valued solutions of the same system, due to 
the invariance under multiplication of either x or z by a cube root of 1. Let 



(1.24a) * := 

(1.24b) ipj(x, z) := tp(x, u> j z), lu := , j = 0,1,2, 

be the Wronskian matrix formed from the three solutions so obtained. Then ^ 
satisfies 

(1.25b) — - 



dz 




which shows that, viewing x as a deformation parameter, the monodromy of the 
rational covariant derivative operator entering in (1.25b) is independent of x. Of 
course, the monodromy in this case is actually trivial, since the matrix Baker func- 
tion is globally defined and single-valued. Moreover since the dependence on 
both the variables x and z is rational, the system could equally well be viewed as 
determining a 1-parametcr family of rational covariant differential operators in the 
x variable, defined by eq. (1.25a), having trivial monodromy, with z viewed as the 
deformation parameter. A dual isomonodromic system may also be defined, just by 
introducing the dual matrix Baker function as the Wronskian matrix determined 
from (^0,^1,^2) through differentioation with respect to the z-variable, 



(1.26a) * := 



(1.26b) ipj(x,z) := iI>(x,uj 3 z), Lu:=e—, j = 0,1,2 



This then satisfies the system obtained from (1.25a,b) by interchanging the x and 
z variables. 

The question naturally arises: is it possible in general to embed bispectral sys- 
tems into isomonodromic ones, such that the matrix Riemann-Hilbert problem used 
in the dressing method construction of solutions to the latter systems reduces to 
the Darboux transformations used in constructing the former? In the following sec- 
tions, this idea will be further developed. A set of criteria is given, within the 
Riemann-Hilbert problem setting, giving rise to isomonodromic systems derived 
from suitably defined "vacuum" solutions. Imposing further restrictions on the 
data entering the Riemann-Hilbert problem, we deduce sufficient conditions for 
the solution to depend rationally on the deformation parameter. Although a set 
of sufficient conditions for these systems to reduce to bispectral ones is not yet 
determined, we are able to illustrate via examples and elementary calculations that 
rank-1 bispectral systems can be recovered from this approach 
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2. ISOMONODROMIC SYSTEMS AND THE RlEMANN— HlLBERT PROBLEM 

2.1. The Dressing Method for Isomonodromic Systems. In the matrix 
Riemann-Hilbert (RH) problem approach to zero-curvature equations depending 
upon a spectral parameter z, one usually starts with a particular r x r matrix 
function ^fo(t,z), viewed as the "vacuum" solution, depending on a number of 
commuting flow parameters t = {ti,t 2 ,- ■ ■) and chooses, along some suitably de- 
fined, closed contour T in the z plane, with interior region T + and exterior T_, a 
smooth, nonsingular r x r matrix-valued function H(z) of the spectral parameter. 
(In the group theoretical formulation, H is interpreted as a loop group element.) 
The Zakharov-Shabat dressing method [NZMP] then consists of solving the associ- 
ated matrix RH problem for the conjugated matrix 

(2.1) x; 1 ^, z) X -{t, z) = * (t, z)H(z)%\t, z), 

where \± are holomorphic functions of z in the regions r± which, in the regular case, 
are nonsingular r x r matrices in their domain of holomorphicity (with a suitable 
normalization condition at some point, usually z = oo, to guarantee uniqueness of 
the solution). The "dressed" matrix Baker function $(x, z) is then obtained by 
applying \± as a gauge transformation to the vacuum solution 

(2.2) * = X +*oH+ = X-^oHZ 1 , 

where the nonsingular r x r matrices H±(z) are t-independent and chosen in any 
way that respects the factorization 

(2.3) H(z) = H+(z)H-(z). 

These just serve to determine a basis for the fundamental solution and may, in 
particular, be chosen either so that eq. (2.3) itself is a solution to the RH problem 
for the fixed element H or, alternatively, that one of the two factors H± is equal to 
H itself, and the other is the identity matrix. 

If the vacuum matrix Baker function satisfies the Zakharov-Shabat (ZS) system 

(2-4) ^ = tf. i0 ( M )* , 

where the C/q's are r x r matrix-valued functions depending rationally on the spec- 
tral parameter z, with a given pole structure and asymptotic limit as z — > oo, 
the dressed Baker function ^(t 7 z) will satisfy a system having the same pole and 
asymptotic structure in z 

(2-5) ^ = ^-(M)*o, 

where the L^-'s are obtained by applying the same gauge transformation to the 
vacuum covariant derivative operators as that appearing in eq. (2.2). It follows that 
the compatibility conditions for (2.4), the zero-curvature conditions, are satisfied 
identically in z. 
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To adapt this procedure to isomonodromic deformation equations, we must as- 
sume that, in addition to the usual vacuum ZS equations (2.4), \P also satisfies a 
compatible differential equation with respect to the spectral parameter 

(2.6) ^=v Q (t,z)^ Q , 

where the rxr matrix- valued function Vq depends polynomially on z. This implies 
that the (generalized) monodromy data of the vacuum rational covariant derivative 
operator 

(2.7) V° z :=-^--V (t,z) 

is invariant under the deformations generated by varying the parameters (ti,t 2 ,- ■ ■)■ 
In order to guarantee that the resulting "dressed" operator 

(2.8) V z :=^-V(t,z), 

obtained through application of the gauge transformation (2.2) to V z , should also 
be rational in the spectral parameter z, further restrictions must be put on the 
choice of H(z). A sufficient set of conditions for H(z) is given in the following. 

Proposition 2.1. If H(z) satisfies a differential equation of the form 

(2.9) =r(z)H(z) + Q(z), 

where r(z) is a rational matrix-valued function and Q(z) is a distribution with 
support at a finite number of points, and the splitting (2.3) is chosen such that 

(2.10) H + (z) := H(z), H_(z)=I, 

the resulting operator V z will depend rationally on z, and its (generalized) mon- 
odromy will be preserved under deformations in the parameters {tj}j=i,2— deter- 
mined compatibly by the equations (2.5). If Q(z) vanishes, then V(t, z) can be 
expressed 

(2.11) V(t,z) = V_(t,z) + V+(t,z), 

where V+ is a polynomial in z, of the same degree as Vq, and V- is rational in z, 
with the same pole support as r{z), and vanishing at z = oo. 

The proof of this proposition is quite elementary and will not be detailed here. 
However, the sense of the derivative appearing in eq. (2.8) needs to be clarified 
since, a priori, H(z) is only defined along the contour T. What is meant here is 
that we assume that H(z) may be extended to a neighborhood of T, though not 
necessarily as a holomorphic function. If it can be extended holomorphically, terms 
like Q(z) are absent. In particular, H(z) could be piecewise constant along T, 
and Q(z) could have its support just on T itself. This case will not be the one of 
interest to us in the following, but it does play an important role in the solution 
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of certain classes of isomonodromic deformation equations through the Riemann- 
Hilbert problem method (cf. [HI]). 

In the following, we shall restrict ourselves to the case when Q(z) = 0, and seek 
to impose suitable restrictions on the vacuum solution ^>q and the element H in 
order that the resulting matrices {Uj(t, z), V(t, z)) appearing in eqs. (2.5) and (2.8) 
depend rationally, not only on the spectral parameter z, but also on the deformation 
parameters {tj}. We begin in the next subsection by giving two examples that lead 
to such rational solutions, and which may be reduced to particular cases of Wilson's 
rank-1 bispectral systems. 

2.2. Two Examples Reducing to Rank-1 Bispectral Systems. 

The following examples show how an appropriately chosen pair (^q, H) leads, 
through solution of the Riemann-Hilbcrt problem, to isomonodromic systems that 
reduce to bispectral ones. We restrict ourselves to 3 x 3 matrices, with \- rational, 
and the pole support of both \~ an d HzH 1 at z = 0. To obtain a correspondence 
between the scalar and matrix systems, we restrict ourselves to the subgroup of the 
full loop group 61(3) consisting of elements satisfying the invariance condition 

(2.12) THiujz)^ 1 = H{z) 

where 



1 











to 











u 2 



(2.13) uj:^e— , T 

with similar rectrictions on Vq and Uq. Define the matrix 
(2.14) 




In both the following examples we choose the vacuum Baker function as 

(2.15) * (x,z):=e xz *, 
so that 

(2.16) U := zE V := xY,. 

From any 3x3 matrix Baker function satisfying the invariance condition (2.12), 
we obtain three scalar Baker functions by summing along the rows; i.e., by taking 
the components of the vector obtained by applying * to the vector with all entries 
equal to 1. 



(2.16) 
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Example 2.2.1. Choose 

(2.17) H(z) : 



as the loop group element. The solution to the Riemann-Hilbert problem is then 
given by 




(2.18) 



1 

i_ 

xz 
2 

x 2 z 2 




which is rational in z. The corresponding isomonodromic system (1.12a,b), (1.13) 
is then defined by 



(2.19) U(x,z)= \ z , V{x,z) 





Summing over the rows of gives 

(2.20) <Pi=e xz , 4> 2 = e xz (l--^J, = e xz (l - ^ 

These are all rank-1 bispectral wave functions which are eigenfunctions, respec- 
tively of the following pairs of operators 



(2.21a) 
(2.21b) 



9 a d 

ox az 



0X Z X z OZ z 



Z 



2 ' 



(2 211 L - — — — — A -— — — — 

3 dx 3 x 2 dx x 3 ' 3 dz 3 z 2 dz z 3 ' 

In terms of Wilson's correspondence, ip\ is just the vacuum Baker function, -02 cor- 
responds to the point in the 1 -particle Calogero-Moser phase space with (X, Z) = 
(0, 0), and tps is just the example of section 1.3, corresponding to the point in the 
2-particle phase space with matrices (X, Z) given in (1.21). Note, however, that 
the correspondence between the isomonodromic and bispectral systems given here 
is not the same as the one given in section (1.3). 

Example 2.2.2. Now choose 
(2.23) H(z) := 
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The solution to the Riemann-Hilbert problem in this case gives 

/ 1 s 

(2.24) X-=\-h 1 

V o -± i 

\ XZ ' 

and the corresponding isomonodromic system (1.12a,b), (1.13) is denned by 

(I z o \ /§ x 

(2.25) [/(z,z) = -I z , V(x,z) = -i x 

\z -i / U -± 

The associated bispectral Baker functions obtained by summing over the rows of ^ 
are 

(2.26) 4>i=e xz , <p 2 = e xz (l-—), ^ = e xz (l- — 

\ xz J \ xz / 

so ipi andV>3 coincide with as ipi an d i>3, respectively, of the previous example, 
while ip2 here is a bispectral eigenfunction of the pair of operators 

(2 27) A- — -— — 

dx 3 x 2 dx ' dz 3 z 2 dz 

This case corresponds to the point in the Calogero-Moser phase space with ma- 
trices 

<**> x-(sj), 

In both these examples, x_ is rational, and its only pole, at z = 0, is coincident 
with that of HzH^ 1 . In the last subsection, some preliminary results will be given 
showing how isomonodromic systems may more generally be obtained, such that the 
dependence on the deformation parameter is also rational, but with an arbitrary 
number of poles in r(z). The structure of the solutions wil be seen to closely 
resemble that of bispectral systems. 

2.3. Rational Isomonodromic Systems. 

We now make some more restrictive assumptions about the solution to the 
Riemann-Hilbert problem, which arc sufficient to assure that the matrices Uj(t,z) 
and V(t, z) are rational functions of the deformation parameters (x — ti, t<i ■ ■ ■ ). 
In the following, we assume that the contour T is sufficiently large that all poles of 
r(z) are contained in the interior region r+, and choose the normalization of X-( z ) 
to be such that 

(2.29) X-(oo) = I. 

We must of course first require that the matrices Uoj(t, z), Vo(t, z) entering in the 
vacuum equations (2.5), (2.6) should themselves depend rationally on the defor- 
mation parameters. Now also assume that the solution to the Riemann-Hilbert 
problem is such \- i s rational in z. It follows that its pole support is contained in 
that of r(z). We now add the additional requirement that the difference 

(2.30) x+(z)^o(z)r(z)^ 1 (z) X + 1 (z) - r(z) 

be holomorphic in T + . This means, for example, when r{z) has only first order 
poles at the points z = Z\, ■ ■ ■ z n , with residue matrices {ri}i—\... n , that the matrix 
X+{zi)^a{zi) must commute with r i7 for each z — Zi. We then have the following 
result 
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Proposition 2.2. Under the above assumptions, (i.e., that X-, Uo_j and Vq are 
all rational in the deformation parameters, and x+ ( z )*&o(z) is such that the ex- 
pression (2.30) is holomorphic in z inside T + ), the matrices Uj(t,z), V(t,z) ob- 
tained through the dressing transformation are also rational in the parameters t = 

(tl,*2"-)- 

Again, the proof is elementary and will not be detailed here. However, to show 
the close similarity with bispectral systems, at least for the case of rank-1 bispectral 
algebras, we will consider the detailed structure of the resulting Baker matrix ^ in 
the case when r(z) has only simple poles and Vo(x, z) is chosen as independent of 
z. (Here x — t\ and the other deformation parameters are chosen to vanish.) This 
means that the vacuum Baker matrix \l/o, as in the scalar case, is a linear exponential 
in z, with matrix exponent zVo(x). (These latter restrictions are imposed only 
for the sake of simplicity and comparison with the bispectral case; the following 
formulae can easily be modified to cover the general case with higher order poles 
and polynomial dependence on z.) For this case, r(z) and X- {x, z) have the general 
form 

n 

(2.31a) rW = £ * 

■ -, z z i 
.7=1 

(2.31b) x _ {XjZ) = i + J2-^, 

where the residue matrices {r^} are constants and the QiS are rational functions 
of x, whose explicit form will be determined under the assumptions of Proposition 
2.2. The relation (2.29) in this case implies the equalities 

(2.32) x+(z l )^o(zr)r l % 1 {z l )x+ 1 (z l ) = r l , i = l,---n. 

It then follows from the dressing method that the polynomial (in z) and rational 
parts V+, V-, respectively, of the matrix V(x, z) entering in Proposition 2.1 are 
just 

(2.33) V+(x,z) = V (x), V-{x,z) = r{z). 

Thus, under the given assumptions, the residue matrices in the polar part of V(x, z) 
are constants, and the polynomial part is the same is for the vacuum solution. The 
dressing transformation gives 

(2.34) (V {x) + r(z)) X -(x, z) = z) + X -(x, z)V (x). 

Substituting the expression (2.31a,b) for X -(x, z) and r(z) and equating the polar 
parts gives the realations 

Qj, j = !•••«• 



(2.35b) 



(2.35b) 



ad 



V 



3 = 1 



*3 
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To express the solution of this linear algebraic system it is convenient to introduce 
the nr x r matrices formed from a column of r x r blocks consisting of the residue 
matrices 



(2.36) 



R:~- 



/ri 



V, 



(Qi 



\Qr 



as well as the nr x nr matrix X, consisting of r x r blocks {Xij}ij—i r .. n given by 

n 

(2.37a) Xii :=-J2-^ 7 

Zi Zj 



(2.37b) 



Zi Zj 



Eq. (2.35a) may then be written more concisely as 



(2.38) 



ady <g>I + X \ Q = -R. 



Solving, and substituting into the expression (2.31b) for X- an d (2.2) for ^ gives 

(2.39) *(x,z)= (l + E T {Z-zl)- 1 {X + a,d Vo ^l)- 1 R S j, 

where E is the nr x r matrix consisting of a column of n blocks, each of which is 
the r x r identity matrix 



(2.40) 



E = 



while Z is the np x np block -diagonal matrix whose r x r diagonal blocks are just 
-zil 

(2.41) Z = -di&g(z 1 ,z 2 ,---z n )®I. 
The commutator of the matrices X and Z is given by 

(2.42) [X, Z] = -Tl + ER T 1 

where TZ is the nr x nr block diagonal matrix with j'th r x r diagonal block equal 
to Rj. 

We see that the rational x-dependence contained in formula (2.39) comes entirely 
from the factor (X+ady C>5l) _1 . Comparison with the formulae in [W] for the Baker 
functions for rank-1 bispectral algebras shows the close resemblance to cq. (2.39), 
while eq. (2.42) is the analogue of the relation (1.5) satisfied by the associated 
Calogero-Moser matrices. 
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Although the above scheme does not guarantee reduction to scalar bispectral 
systems, the examples of the previous subsection show that they do so at least in 
particular cases. Also, although the conditions of Proposition (2,2) can be used 
constructively when looking for isomonodromic systems involving rational depen- 
dence on the deformation parameters, it would be more satisfying to find sufficient 
conditions for rationality, and for bispectrality, purely in terms of the input data 
consisting of the vacuum solution "to and the group element H . 



Acknowledgements. The author would like to thank A. Kasman and A. Orlov for 
helpful discussions. 
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